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ABSTRACT 

We consider the dependence of ultra-high energy cosmic ray (UHECR) flux predic- 
tions on the choice of galax y catalogue. We demonstrate that model predictions by 
iKoers fc Tinvakovl ()2009bD . based on the so-called KKKST catalogue, are in good 
agreement with predictions based on the XSCz catalogue, a recently compiled cat- 
alogue that contains spectrosco pic redshifts f or a large fraction of galaxies. This 
agreement refutes the claim by iKashtil (120091 ) that the KKKST catalogue is not 
suited for studies of UHECR anisotropy due to its dependence on photometric 
redshift estimates. In order to quantify the effect of galaxy catalogues on flux pre- 
dictions, we develop a measure of anisotropics associated with model flux maps. 
This measure offers a general criterion to study the effect of model parameters 
and assumptions on the predicted strength of UHECR anisotropies. 
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1 INTRODUCTION 

Galaxy catalogues are an indispensable tool to model the 
structure of the local Universe, as required for studies of 
ultra-high energy cosmic ray (UHECR) anisotropy. Many 
such studies (e.g. Cuoco et al.l 20061: Kashti fc WaxmanI 



l2008l : iTakami et al.1 l2009l: iKa shti' '2009() have used the 
PSCz catalogue l|Saunders et al. 2000) for this purpose. 
This catalogue, however, has its drawbacks: it suffers from 
incomplete sky coverage, it m ay underestimate g alaxy 
counts in high-density regions (|Huchra et al.|[2007h . and 
i t has limited stat i stics. The 2MASS galaxy catalogue 
(IJarrett et al ] |200d : lJarrettll2004l : [skrutskie. M. F. et al.l 
I2OO6I ) . offering complete sky coverage (except for the 
galactic plane) and excellent statistics, improves on these 
issues. However, this catalogue does not contain spectro- 
scopic redshifts so that redshifts have to be estimated by 
photometry, i.e. from the observed brightness under the 
assumption of a standard intrinsic brightness in a spe- 
cific band. Efforts to obtain spectroscopic redshifts for 
2MASS galaxies are ongoing. The resulting catalogue, 
termed XSCz, is presently being compiledQ 

Mot ivated by the imperfec tion of existing galaxy cat- 
alogues, IXalashev et al.l (|2008l ) have compiled a hybrid 
catalogue that contains galaxies with spec troscopic red- 
shifts from the (HYPER)LEDA database l|Paturel et al.l 
I2OO3I ) within 30 Mpc and galaxies from the 2MASS Ex- 



tended Source catalogue (|jarrett et al.l2000l ) at larger dis- 
tances, for which photometric redshift estimates are used 
(the two contributions are of course normalized appropri- 
ately) . This catalogu e, which we have termed th e KKKST 
catalogue, is used in lKoers fc Tinvakovl l|2009bl ; hereafter 
paper I) 



The adequacy of the KKK ST cata l ogue for UHECR 
flux modelling is questioned bv iKashtil l|2009l ) because of 
its dependence on photometric redshift estimates. The er- 
rors in these estimates give rise to errors in flux estimates, 
which may distort model flux maps and may lead to inac- 
curate model predictions. The uncertaint y in pho t omet - 
ric redshift estimates is indeed fairly large: Ijarretti (|2004l ) 
quotes an error of 10% — 20% for mos t normal galax- 
ies in the 2MASS survey (|Kashtil (|2009l ) claims a ~30% 
systematic uncertainty). In paper I, we deemed these un- 
certainties acceptable for two reasons. First, large-scale 
anisotropies in flux maps arise as a collective effect of 
many sources. Adding flux contributions of many indi- 
vidual galaxies reduces the relative strength of fluctua- 
tions (for the KKKST catalogue with a smearing angle of 
&s = 6°, as adopted in paper I, a model flux is composed 
of individual contributions of ~600 galaxies). Second, the 
remaining uncertainties only affect the model flux beyond 
30 Mpc, while the imprint of local structure is strongest 
at close distances. 



* E-mail: hkoers@ulb.ac.be 

f E-mail: Petr.TiniakoviOlulb. ac.be 

^ See http://web.ipac.caltech.edu/staff/jarrett/XSCz/. 



In this paper we compare model fluxes based on the 
KKKST catalogue to model fluxes based on a preliminary 
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version of the XSCz catalogue0 This aUows us to deter- 
mine to which extent the flux predictions in paper I are 
contaminated by the inaccuracies in the catalogue used. 
The XSCz catalogue is very well suited to model the mat- 
ter distribution in the Universe because of its complete- 
ness, large statistics, and spectroscopic redshift measure- 
ments. We therefore consider the resulting model fluxes 
as a benchmark against which UHECR flux predictions 
derived with the KKKST catalogue (as well as with the 
PSCz catalogue) can be cross checked. As we will demon- 
strate, flux predictions based on the KKKST catalogue 
(as well as on the PSCz) are in good agreement with the 
XSCz prediction. This provides an a posteriori verifica- 
tion of the accuracy of our models in paper I. 

We also investigate, in a general setup, the rela- 
tionship between galaxy catalogues and the predicted 
strength of UHECR anisotropics. For this purpose we 
develop a measure that quantifies the relation between 
model flux maps, which represent the model predictions 
for a given galaxy catalogue, and the test statistic D (in- 
troduced in paper I), which is a measure of the strength of 
UHECR anisotropics. This measure can be used to assess 
the predicted strength of anisotropics from the "contrast" 
that is exhibited in model flux maps and is independent of 
event number. Because of its generality, this measure may 
also be useful to study the effect of other model param- 
eters and assumptions (e.g., threshold energy, deflection 
angle, or injection spectrum) on the predicted strength of 
UHECR anisotropics. 

This paper is organized as follows. In section [5] we 
compare model predictions from the KKKST and PSCz 
catalogues to predictions based on the XSCz catalogue. 
Section O concerns the relationship between galaxy cata- 
logues and the predicted strength of UHECR anisotropics. 
We summarize our flndings in section |4l 



2 ACCURACY OF MODEL FLUX MAPS 

We begin our analysis by considering the accu- 
racy of photometric redshift determinations in the 
2MASS catalogue. We compute the relative difference 
5z/z :— (2:photo — Zspcc)/2:spGc, where ziphoto denotes a pho- 
tometric redshift estimate and Zspec denotes a spectro- 
scopic redshift, for all galaxies in the XSCz catalogue that 
have spectroscopic redshifts. Here, and throughout the 
paper, we only consider galaxies with /sTs-magnitude be- 
low 12.50 in the XSCz catalogue. Spectroscopic redshifts 
are available for ^^70% of these relatively bright galaxies. 
The distribution of Sz/z is shown in Fig. [T] The inac- 
curacy in photometric redshift estimates is indeed fairly 
large: we find that the average value of \5z/z\, the abso- 
lute relative difference in redshift determinations, is 0.36. 

The UHECR luminosity at the source being un- 
known, model fiux predictions are conveniently formu- 
lated in terms of the normalized flux 

^ We thank Tom Jarrett for providing us with a preliminary 
version of the XSCz catalogue. 
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Figure 1. Distribution of the relative difference 5z/z between 
photometric and spectroscopic redshift estimates for individual 
galaxies in the XSCz catalogue. Inset: cumulative distribution 
of \5z/z\. In producing this figure we only used galaxies with 
D > 5 Mpc. 
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Figure 2. Distribution of Af for individual galaxies in the 
XSCz catalogue. In producing this plot, we approximated 
fluxes as oc and we only considered galaxies with D > 5 
Mpc. 

where (p denotes a single-source flux and {</>) is the average 
flux of all galaxies. (Note that the overall flux normaliza- 
tion is determined by observations). To get a rough esti- 
mate of the inaccuracies in flux predictions at the level 
of individual sources, we approximate the flux associated 
with every galaxy as oc D~^, D being the source dis- 
tance (we take D on z here). Here, and throughout this 
paper, we assume all sources have the same intrinsic lumi- 
nosity. We then compute the quantity A f = /photo — /spec, 
where /photo is the normalized flux estimate based on the 
photometric redshift estimate and /spec is the normalized 
flux derived from the spectroscopic redshift. The distri- 
bution of Af is shown in Fig. [2] 

We now consider model UHECR flux maps, which 
are generated using the methods described in paper 
I. Throughout this work we assume a proton compo- 
sition and a power-law injection spectrum with spec- 
tral index p = 2.2 extending to very high energies. We 
adopt a ACDM concordance model with Hubble constant 
72 km s~^ Mpc~^, and cosmological density parameters 
Qm ~ 0.27 and Qa ~ 0.73. Following the setup in paper 
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Figure 3. UHECR model flux maps for Eq = 60 EeV, = 6°, 
and uniform exposure. The three panels show model predic- 
tions based on different galaxy catalogues: the KKKST cata- 
logue (top panel), the XSCz catalogue (middle) and the PSCz 
catalogue (bottom) . The gray bands are chosen such that each 
contains 1/5 of the model flux, with darker gray indicating 
larger flux. 



I, we adopt a smearing allele ffa = 6°. UsinR the method 
described in lKoers fc Tinvakovl lj2009ah . we have verified 
that all catalogues contain sufficient galaxies to provide 
a good statistical description on this angular scale. We 
remove the region |6| < 15° from our analysis because of 
incompleteness near the galactic plane. (We choose the 
same region for every catalogue for comparison; for the 
XSCz skymap this choice is overly conservative). Fig. [3] 
shows the flux maps for a threshold energy Eo = 60 
EeV. In Fig. |4] we show the distribution of flux differ- 
ences A/ — /alt — ./xscz, where / — 4>/{4>) denotes the 
normalized flux and the subscript refers to the catalogue 
that was used in the modelling ( "alt" standing for either 
KKKST or PSCz). Here (f) represents the integral UHECR 
flux in a given direction (as represented in Fig. |3]) and (0) 
is the average flux on the sphere. 

Figures [3] and |4] demonstrate that model predictions 
based on both the KKKST and PSCz catalogues are in 
good agreement with those from the XSCz catalogue. The 
average value of |A/| is 0.17 for the KKKST catalogue 
and 0.19 for the PSCz. At lower threshold energies the 
differences become smaller, while at higher energies they 
are somewhat larger. At _Eo = 100 EeV, for example, the 
average |A/| is 0.39 for the KKKST catalogue and 0.32 
for the PSCz. 



Sampling a flux map uniformly over the sky, one ob- 
tains a flux distribution associated with the map. This 
distribution encodes the relevant intrinsic properties of 
the flux map, i.e. properties relating to the strength of 
over- and underdense regions but not to their position on 
the sky. Postponing a more thorough discussion of the flux 
distribution to the next section, we point out here that 
the flux distribution corresponding to the KKKST flux 
map is broader than the flux distribution obtained with 
the PSCz catalogue (the XSCz is in between). This is re- 
flected in the fact that the band of highest flux (darkest 
gray) in Fig. [3] occupies the smallest area on the sphere: 
7.5% compared to 8.2% for the XSCz flux map and 10.4% 
for the PSCz flux map. For an isotropic flux map this 
number would be 1/5 = 20%. We see that the KKKST 
map deviates most strongly from an isotropic flux map. As 
a consequence, model predictions based on this catalogue 
will exhibit the strongest departures from isotropy. We 
will use the term "contrast" to refer to the width of the 
flux distribution, strong contrast corresponding to a wide 
flux distribution and a strong deviation from isotropy. 

How does the choice of galaxy catalogue affect the 
statistical tests proposed in paper I? We consider here 
the D-test, to which we will refer as the "flux sam- 
pling method". The test deflnes a test statistic D as the 
Kolmogorov-Smirnov (KS) distance Dks between the cu- 
mulative distribution of flux values sampled by a set of 
UHECR events and a reference distribution that corre- 
sponds to the model that is tested (see paper I for de- 
tails and discussion). Here and in the following we con- 
sider two models: the "Isotropy" model (denoted I), which 
states that UHECR events are distributed isotropically 
(we do not consider experimental exposure here), and 
the "Structure" model (denoted S), which states that 
UHECR sources trace the distribution of matter in the 
Universe. As a case study, we show in Fig. [5] the distri- 
bution of the test statistic D for 20 UHECR events with 
energies in excess of 60 EeV. Here the reference distribu- 
tion corresponds to an isotropic flux, so that the curve la- 
beled "Isotropy" follows the universal KS self-correlation 
distribution (it is thus the same for all catalogues). The 
"Structure" distribution obtained from the PSCz cata- 
logue is clearly closer to the "Isotropy" distribution than 
the "Structure" distribution obtained from the KKKST 
catalogue (the distribution from the XSCz being interme- 
diary). This is reflected in the statistical power estimates: 
For a significance a = 0.05, we find P = 0.54 for the 
KKKST catalogue, P — 0.46 for the XSCz catalogue, and 
P — 0.37 for the PSCz catalogue, where P denotes the 
power to reject I when S is true. The difference in statisti- 
cal powers can be related to the contrast in the model fiux 
maps discussed in the previous paragraph. This relation 
will be explored in the next section. 

We have recomputed power estimates and p-values 
for the D-test reported in paper I with the XSCz cata- 
logue instead of the KKKST catalogue. At energies below 
60 EeV we find that the power estimates change by only 
a few percent. At 100 EeV the XSCz catalogue leads to 
power estimates smaller by ~20%, which is still accept- 
able in the light of other uncertainties. For the data ob- 
tained by the Pierre Auger Observatory and the AGASA 
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Figure 4. Distribution of Af for the model fluxes shown in 
Fig. [3] 
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Figure 6. The cumulative distribution Cisoif) for the KKKST, 
PSCz, and XSCz catalogues. This figure applies to threshold 
energy Eg = 60 EeV and smearing angle 9s = 6° . 
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Figure 5. Distribution of test statistic D for an isotropic ref- 
erence distribution and test UHECR event sets following the 
distribution of matter as modelled from different catalogues. 
This figure applies to 20 events with threshold energy 60 EeV, 
and the distributions are obtained by sampling over 10* inde- 
pendent Monte Carlo realizations. 

experiment, we find no significant difference in p- values 
for the different catalogues. 



3 GALAXY CATALOGUES AND 
PREDICTED STRENGTH OF 
ANISOTROPIES 

In this section we investigate how the choice of cata- 
logue affects the predicted strength of large-scale UHECR 
anisotropies. In order to quantify this strength we use the 
fiux sampling test to measure deviations from isotropy. 
Although other tests will measure anisotropies in a dif- 
ferent manner, a galaxy catalogue that yields a value for 
the D test statistic close to the isotropic prediction will in 
general also yield outcomes close to the isotropic predic- 
tion in other tests (and similarly for predictions far away 
from the isotropic one). 

The outcome of the fiux sampling test depends on 
the number of events as well as on the distribution of 
model fiuxes over the sky, i.e. the fiux map. The first pa- 
rameter is obviously independent of the choice of galaxy 



catalogue. All intrinsic properties of a galaxy catalogue 
that may affect the strength of anisotropies are contained 
in the fiux map. Here and in the following, the expression 
"flux map" is used to refer to the flux map derived under 
S (i.e., assuming UHECR sources trace the distribution 
of matter in the Universe) . The flux map also depends on 
threshold energy, UHECR injection spectrum, and aver- 
age deflection angle. 

The connection between flux maps and the strength 
of UHECR anisotropies can be illustrated by two limit- 
ing cases. First, consider the case that S would predict 
a uniform fiux, so that the fiux distribution would be a 
delta- function (minimal contrast). In that case UHECR 
events sample the sky uniformly under both S and I and 
model predictions become identical. In the opposite case 
of a very wide fiux distribution (strong contrast), events 
will have a strong tendency to cluster in high-fiux regions 
under S, exhibiting strong anisotropy. The connection be- 
tween fiux distributions and the strength of anisotropies 
can be made explicit using the D-test. As shown in ap- 
pendix [X] under S the test statistic D approaches the 
following limiting value as the number of events goes to 
infinity (recall that it approaches in the same limit un- 
der I): 

f dfCi^^if), (2) 
Jo 

where Ciso(/) denotes the normalized cumulative fiux dis- 
tribution under S and we recall that / = (j> /{())), (j) being 
the integral UHECR fiux and {(j)) the average value on the 
sphere. Equation ((2)1 identifies the intrinsic property of a 
fiux map that is important in determining the strength 
of UHECR anisotropies. It gives a precise meaning to the 
concept of "contrast" that was introduced in the previous 
section. It is related to the width of the fiux distribution: 
A narrow fiux distribution corresponds to a steeply rising 
cumulative distribution and, via equation ([2]), a value of 
D close to 0. 

The distribution Ciso is shown in Fig. for the three 
fiux maps presented in Fig. [S] Comparing the surfaces 
under the curves for < / < 1, it is clear that the PSCz 
catalogue yields the smallest asymptotic value of D, and 
the KKKST catalogue the largest. The larger contrast 
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Figure 7. Asymptotic value of D for A^cv oo for the 
KKKST, XSCz, and PSCz catalogues as a function of Eq. 



in the KKKST catalogue implies stronger signatures of 
anisotropy. This explains the ordering of statistical powers 
that was found in the previous section. 

In Fig. [7] we show the asymptotic value of D from 
eq. ([U for different catalogues and threshold energies. We 
observe that the KKKST catalogue systematically pre- 
dicts somewhat stronger anisotropics compared to the 
PSCz catalogue, in keeping with the results presented 
above and with the discussion in paper 10 The XSCz 
prediction lies in between. For threshold energies Eq < 60 
EeV, the XSCz curve is virtually identical to the KKKST 
one, while it approaches the PSCz curve at higher ener- 
gies. The differences between the three curves are how- 
ever small. In particular, they are smaller than the un- 
certainty induced by systematic errors in energy determi- 
nation. (One may check that the curves corresponding to 
the KKKST and PSCz catalogues are within a 20% shift 
in energy applied to the XSCz prediction). 

In reality the number of UHECR events is of course 
finite, and often not very large. The discriminatory power 
then depends not only on the flux map but also on event 
number. In this case, equation provides a figure-of- 
merit that governs the asymptotic value of D. This is 
illustrated in Fig. [51 where we show the average value 
of D and the l-cr bands as a function of event number. 
When A^cv becomes very large, D approaches the asymp- 
totic value computed with eq. ©: 0.22 for the KKKST 
catalogue and 0,17 for the PSCz catalogue. 



4 SUMMARY 

In this paper we have investigated how the choice of 
galaxy catalogue affects UHECR model fluxes in a sce- 
nario where UHECR sources trace the distribution of 



^ Unfortunately, the power estimates for the PSCz catalogue 
as reported in paper I are inaccurate due to an erroneous appli- 
cation of the PSCz selection function in our numerical routines. 
In table 3 the power estimates for the PSCz catalogue should 
read 0.37, 0.28, and 0.72 for scenarios I, II, and III, respec- 
tively. With these changes, the statistical powers for the PSCz 
catalogue remain smaller than those for the KKKST catalogue 
so that the conclusions remain qualitatively unchanged. 



Figure 8. Average value of D (thick lines) and 1-sigma con- 
tours (thin) as a function of event number N^v for the KKKST 
and PSCz catalogues and for = 60 EeV. The full distribu- 
tion for Ncv = 20 is shown in Fig. [5] 



matter in the Universe. The differences between the 
three catalogues consi dered here, the KKKST catalogue 
(|Kalashev et all 120081) . the XSCz cata logufl and the 
PSCz catalogue (jSaunders et ahl |2000| ) are reasonably 
small. This is reassuring because all catalogues are sup- 
posed to be sampling the same underlying density field 
(barring biases induced by selection effects). 

In section [2] we have shown that there is good agree- 
ment between model flux maps constructed with the 
KKKST and XSCz catalogues. The former wa s used by us 
earlier in paper I (|Koers fc Tinvakovll2009bl ) The agree- 
ment is especially good at energies below 60 EeV, the 
regime where we have confronted models with experimen- 
tal data in p aper I. This c omparison refutes the recent 
statement by iKashtil (|2009l ) that the KKKST catalogue 
is not suited for studies of UHECR anisotropy because of 
its dependence on photometric redshift estimates. 

We have investigated the relation between the pre- 
dicted strength of large-scale UHECR anisotropics and 
model flux maps from a general point of view in section 
[21 The intrinsic properties of an UHECR flux map, i.e. 
those properties relating to the strength of over- and un- 
derdense regions but not to their position on the sky, are 
contained in the flux distribution. Equation ^ demon- 
strates how this distribution can be used to determine the 
value of the D test statistic in the limit of infinite events. 
This asymptotic value provides a measure of the expected 
anisotropy in UHECR arrival directions for sources trac- 
ing the distribution of matter in the Universe. We have 
compared these values for the KKKST, XSCz, and PSCz 
catalogues as a function of energy (see Fig. (Tjl . The com- 
parison shows that the KKKST catalogue typically yields 
stronger anisotropics than the PSCz catalogue, the XSCz 
catalogue being in between. The difference is however 
small in the light of the uncertainties induced by system- 
atical errors in UHECR energy determination. 



See http: //web. ipac , caltech.edu/staff/jarrett/XSCz/ 
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APPENDIX A: KOLMOGOROV-SMIRNOV 
DISTANCE IN THE LIMIT OF INFINITE 
STATISTICS 

In this appendix we derive equation . We choose a gen- 
eral setup because the result is not limited to UHECR 
anisotropy tests, but is applicable to any experimental 
test that can be appropriately formulated. 

Consider an experiment that records events which 
are characterized by a set of quantities x (these quantities 
may be observables or derived quantities). The space of all 
possible experimental results is denoted as V. Now con- 
sider two models, termed model A and model B. Model A 
asserts that events sample V uniformly, i.e. that the prob- 
ability of registering quantities in the range a; ... a; -I- da; is 
independent of x (as long as x £ V). Within model B, on 
the other hand, the probability that an event has quan- 
tities in the range a; . . . a; -I- da; is proportional to f{x)dx, 
where / is a non-trivial, known function. This function 
defines a map a; i— > /, which associates a real number / 
to every event. Note that we leave the normalization of / 
arbitrary. 

Our aim is now to differentiate between models A 
and B. We will do this by considering the distribution of 



/ over an observation, i.e. a series of registered events. We 
consider the limit of infinite statistics. The distribution 
of / under model A, denoted as I>a(/), represents the 
distribution of / over V. We assume this distribution is 
normalized, i.e. 



d/©A(/) = l. 



(Al) 



The distribution of / under model B is denoted as Ob (/) . 
By construction. 



Ob(/)=c/Oa(/), 
where 

c = I 1^ dffVAif) 



(A2) 



(A3) 



is a normalization constant to ensure that 7?b(/) is also 
normalized. Note that 1/c coincides with /a, the average 
value of / in model A. We define cumulative distribution 
functions as 



C(/) 



/ 



df'Vif) 



(A4) 



The Kolmogorov-Smirnov distance Dks is a measure of 
the difference between Ca and Cb. As we will show, it 
can be expressed in terms of the cumulative distribution 
function Ca(/) alone. First we recall that, by definition, 
Dks = max |L»(/)|, where D{f) = Ca(/) -Cb(/) denotes 
the difference between the two cumulative distribution 
functions. This difference can be expressed as follows: 



D{f)= / d/'©A(/')(l-c/') . 



(A5) 



Since ©a, c and /' are strictly positive, the integral is 
maximal when / = 1/c = /a . Hence 



1 



/a 



df'CAif) 



(A6) 



where the last equality has been obtained by integration 
by parts. Normalizing the function / such that /a = 1, 
the r.h.s. of eq. (|A6|I reduces to the integral of the cumu- 
lative distribution function Ca up to / = 1. 

Two remarks are in order. First, equation (|A6|) is a 
general result that is applicable to any statistical test for 
a model that can be formulated in terms of a probability 
function f{x), where x is distributed uniformly in V under 
the alternative model. The result applies to the UHECR 
anisotropy tests considered in this work by associating 
X = {l,amb), where I and b denote galactic coordinates, 
and / = (f>/{(j>}, where (p denotes the model fiux when 
UHECR sources trace the distribution of matter, and {(p) 
is the average value on the sphere. Note that, by defi- 
nition, /a = 1 in this case. Second, equation HA6|I deter- 
mines -Dks in the limit of an infinite number of events. For 
finite event numbers, it is still a useful figure of merit be- 
cause it determines the asymptotic behaviour (cf. Fig. [8)l . 
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